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The study of highly photo-excited matter at solid state density is an emerging field of research, which is 
benefitting the development of free-electron-laser (FEL) technology. We report an extreme ultraviolet 
(XUV) reflectivity experiment from a titanium (Ti) sample irradiated with ultrafast seeded FEL pulses at 
variable incident photon fluence and frequency. Using a Drude formahsm we relate the observed increase in 
reflectivity as a function of the excitation fluence to an increase in the plasma frequency, which allows us to 
estimate the free electron density in the excited sample. The extreme simplicity of the experimental setup 
makes the present approach potentially a valuable complementary tool to determine the average ionization 
state of the excited sample, information of primary relevance for understanding the physics of matter under 
extreme conditions. 

The study of matter under extreme thermodynamic conditions is a longstanding and fruitful field of research 
that has generated a host of applications in fundamental physics and applied cutting-edge technology''^. In 
the last few decades increasing efforts have been devoted to studying new states of matter such as strongly 
coupled plasma and warm dense matter (WDM)''. The latter is characterised by solid-like density and tempera- 
tures of a few eV, conditions actually found in astrophysical contexts* or in devices for plasma creation and 
inertial confinement fusion^. The advent of ultra-fast/ultra- intense lasers has opened up the possibility to generate 
and probe such extreme states of matter in transient conditions^ '. However, the capability of optical lasers to 
probe the excited sample is hampered by the "screening effect" of the excited plasma, which drastically limits the 
penetration length for photons with frequency (ffl) lower than the plasma frequency: 



(1) 



where JV^, e, m,, and £o are the free electron density, the elementary charge, the effective electron mass and the 
dielectric constant, respectively^''. 

In this context, the development of Free Electron Lasers (FELs) operating in the extreme ultraviolet (XUV) and 
x-ray spectral ranges represent a landmark result"*, since the high energy photons emitted by these sources can 
penetrate matter at solid densities ensuring uniform excitation of the sample. FEL-based experiments have 
provided valuable new insights into this field, for example saturable absorption by core electrons at XUV/x- 
ray wavelengths" or depression of the ionization potential in solid-density plasma'-". Other FEL-based experiments 
aimed at the determination of the relevant parameters of these highly excited states of matter have been proposed, 
such as x-ray induced changes in optical reflectivity'" the parallel pyrometrical characterization of heat waves 
for the determination of the ion temperature (T,)'^ or Thompson scattering for the determination of both free 
electron density (N^) and temperature (T^Y*. The latter approach is based on the measurement of the plasmon 
lineshape in the inelastic scattering spectra of the excited sample. These highly informative experiments rely on 
the development of sophisticated XUV/x-ray spectrometers and are limited by coherence issues'^ and, in high-Z 
samples, by the low scattering cross-sections'''. 

Information on the plasma resonance (and thus on the related sample parameters) can be also gained by 
simpler reflectivity measurements""*. For instance, a shift in cUp upon photo-excitation, observable in Thompson 
scattering experiments, can be ascribed, through Eq. (1), to a change in Ng. However, the same shift in ojp also 
leads to a change in the XUV reflectivity. This is illustrated in Fig. 1, which reports the XUV reflectivity (R) of a 
metal at given values of co for a typical range in ojp and calculated using the Drude model (here y is the inverse 
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Figure 1 | ?j»p-dependenceofRforfic;= 18.9 and 20 eV (black/red lines) 
and hy = 1, 1.5 and 2 eV (higher to lower curves), calculated with Eqs.(3) 
and (4); the upper horizontal scale is the corresponding range in JV^ 
calculated through Eq. (1). The vertical blue segment indicates the 
expected value of ScOp = ?!Ojp(0) = 17.7 eV for the unexcited Ti sample"'^*, 
while the hatched area sketches the fitUp-range reached in the excited 
sample in the exploited f-range. 

lifetime of the plasma resonance""^', see Methods for further details). 
If CO is tuned just above lOp then the dependence of R on cOp is more 
significant than that on In this case any variation in R can be 
essentially attributed to a variation in cOp and thus in Ng. Following 
this idea we performed an XUV reflectivity experiment on a Ti sam- 
ple as a function of both photon fluence (F) and cu. The experiment 
was carried out at the TIMEX end-station of the Elastic and Inelastic 
Scattering (EIS) beamline"'^'"'\ exploiting the broad co-tunabUity in 
the low frequency side of the XUV spectrum provided by the FERMI 
EEL source""^''; experimental details are reported in the Methods 
section. 

Results 

The dependence of i? on P is shown in Fig. 2 for three representative 
values of hoj and clearly show an increase of reflectivity up to about 
40-50% for F-values as large as 20-25 J/cm"'. No appreciable changes 
in _R has been observed for F < 0.02 J/cm^ (see the Methods section). 
Assuming that the increase in (and thus that in Wp) is proportional 
to the energy density (£) deposited into the sample close to the metal 
surface, it is possible to calculate the _F-dependence of _R at a given 
value of CO. For this purpose we used the Drude model to account for 
metal reflectivity, in which we consider the time and transverse spa- 
tial profiles of the FEL pulse and the fractional increase in N,, per unit 
E, i.e.: 

ffl; = co;(0)(l+A£), (2) 

where co^(O) = 17.7 eV^''^** is the plasma frequency of the unper- 
turbed sample (more details of the model used for the calculations 
are reported in the Methods section) and A is an empirical constant. 
Curves reported in Fig. 2 were obtained using A = 7.5 • 10'^ ra^/]; 
higher/lower dashed lines correspond to A-value 25% larger/smaller 
than 7.5 • 10'^ m'/J and can be regarded as an estimate of the con- 
fidence interval (see Supplemental Information), which roughly 
matches the deviation of the data from the calculated curves. 

Once the parameter A is determined, it is straightforward to com- 
pute the F-dependencies of cOp and N^, and these indicate an increase 
in tOp and N,, (averaged over both the time and transverse spatial 
profiles of the FEL pulse) of up to =4% and 8%, respectively, at the 
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Figure 2 | F-dependence of R{F)/R{0) for some selected fisv-values, 
indicated in the individual panels. Each data corresponds to a single shot 
measurement taken in a fresh portion of the sample. Red full lines are the 
R(F)/R{Q) trends calculated through Eq. (6); dashed Unes are the estimate 
of the confidence interval. 

highest F-values. Fig. 3 reports the variation with F of the mean ion 
charge: Z(F) = Z(0)]V,(F)/N,(0), where Z(0) = 4"-"'. Such a variation 
is most likely due to the ionization of the shallower (3p) core elec- 
trons, that cannot efficiently recombine within the FEL pulse dura- 
tion, thus leading to a net increase in the average density of free 
electrons. We note that the F-values reached in the present experi- 
ment compares well with the ionization energy (F* = 99 eV/ 
atom^'''") for the Ti"*^ Ti^* process. If we assume that the latter 
is the leading ionization process, then the data shown in Fig. 3 indi- 
cates that at the highest probed Fs and within the —45 fs FEL pulse 
duration =40% of Ti atoms are ionized. This is in qualitative agree- 
ment with the results reported in Ref. 8 for Al, where the onset of a 
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Figure 3 | Full lines are the F-dependencies of Z(F) at the employed ho> 
values; dashed Unes are the confidence range. Inset compares the low F 
estimate of F, with data from Ref. 8 (circles). 
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Figure 4 | Sketcli of the experimental setup (side view - not to scale): GA, GM, SM and PD are the gas attenuator cell, the gas monitor, the spherical 
focusing mirror and the photodiode, respectively. The FEL beam is shown in green. The relevant distances and angles are indicated. 



FEL-induced saturable absorption is observed in a similar F-range. 
This corresponds to the situation where almost all atoms are 
ionized within the FEL pulse duration. Furthermore, by mak- 
ing the crude assumption that in a moderately excited sample 
£ = (3/2)fcB T,(F)Z(F) + (Z(F) -Z(0))£+' one can derive T, from 
Z. The values estimated in this way are shown in the inset of Fig. 3 
and compared with those reported in ref 8 for Al. The latter data were 
obtained by interpreting the soft x-ray emission spectra within a free 
electron gas model (like the one used in the present work), so that we 
can speculatively infer a common low-F behavior of highly photo- 
excited metals in the sub-0.1 ps timescale, which is mainly deter- 
mined by the role played by free electrons. 

Discussion 

We have reported the experimental observation of a large XUV 
reflectivity enhancement in titanium upon irradiation with ultrafast 
high-fluence FEL-pulses, whose photon frequency was tuned above 
the plasma frequency of the sample. The FEL irradiation drives the 
titanium sample into a transient (short-lived) excited state that exhi- 
bits electronic temperatures typical of dense plasmas but still retains 
a cold atomic lattice. Reflectivity data were interpreted through a 
simple analytical model, based on the free electron gas approach, 
which employs a single empirical parameter to account for the 
increase in the plasma frequency within the FEL pulse duration. 
Such a simplified data analysis is able to capture the main physical 
process at the base of the observed phenomenology and allows us to 
determine the mean free electron density of the excited sample. 

The information that can be gained from the present experiment is 
of the highest relevance for the study of highly excited states of 
matter. Analogous information can be deduced from the lineshape 
analysis of Thompson scattering spectra. However, the present 
approach has the advantage of the extreme simplification of the 
experimental setup and, additionally, reflectivity measurements do 
not suffer from cross-section or coherence issues. Substantial 
improvements in the experimental accuracy and data analysis can 
be achieved with moderate efforts and will likely open up the pos- 
sibility to obtain additional information from this class of experi- 
ments, for example possible modifications of the electronic band 
structure. Also, the capability of FERMI to operate in the near future 
at photon energies down to —12 eV will allow one to probe the 
condition where the metal reflectivity is almost insensitive to plasma 
frequency variations but strongly depends on the plasmon lifetime. 
Since the latter is mostly determined by electron-electron and elec- 
tron-phonon interactions, information on these dynamics could be 
gained as well. In summary, the possibility to map out FEL-induced 
reflectivity variations in a broad range across the low frequency side 
of the XUV spectrum can be profitably used as a complementary tool 
with respect to Thompson scattering methods by the growing com- 
munity involved in the study of matter under extreme conditions. 



Methods 

XUV reflectivity measurements. A sketch of the experimental setup is shown in 
Fig. 4. A Ti mirror (substrate: Si, roughness ~1 nm I^S, thickness 100 nm, 
passivated with 3 nm Ti02) was loaded into the 5-axis manipulator (x,y,z,pitch,roll) 
of the EIS-TIMEX end-station. The FEL beam was focused onto the sample by a 
spherical platinum-coated silicon mirror (400 mm radius of curvature, 0.2 nm 
roughness RMS) placed close to normal incidence (angle of incidence 3'^), which is 
able to provide a focal spot of area o"^ ~ 100 /.im^. The maximum energy/pulse 
delivered by the FEL source was in the —300 //J level. The actual energy/pulse at the 
input of the photon transport system of the EIS-TIJVIEX beamline {I pel) was varied in 
the 1-300 fi] range by a gas attenuator and monitored on a shot-to-shot base by a gas 
monitor^\ Each shot was tagged by a unique label (bunch number) to one-to-one 
correlate all the acquired data. F was determined from Ifj^i taking into account cr^ as 
well as the transmission of the Ti02 layer, the beamline transmission and the mirror 
reflectivity that, in the used oj-range, are in the —73%, 52-56% and 15-20% ranges, 
respectively. Appreciable damage of the sample surface was observed after irradiation 
with a single FEL pulse with f > 0.1 J/cm^. The pulse energy reflected by the sample 
(lit), at 6" angle of incidence, was collected at ho> = 18.9, 19.1, 19.4, 19.7, 20.0 and 
20.3 eV by a Silicon photodiode (UVG20S, IRD inc) coupled with a 0.5 mm thick 
YAG fluorescence screen having a 100 nm aluminum coating on the FEL side. In 
order to calibrate the high-f experiment we collected the ratio Ir/Ifei —Ir/ Ifel in the 
1-300 /fj /^■j./^-range after having moved the sample 3 mm away from the focal plane 
(here F < 0.02 ]/cm^). The results are shown in Fig. 5 and display a constant ratio, 
indicating no appreciable changes in the reflectivity in this low-F range. The I^/lfEL 
value was used to determine the relative reflectivity variation in the high-F 
experiment, i.e.: R{F)/R(0) = (/r//fel)/(Jr/4i)- 

Data analysis. The reflection coefficient at near normal incidence can be expressed in 
terms of the complex dielectric function (c) as^: 



(3) 



which is a sufficiently accurate approximation for sample thickness much larger than 
the penetration depth of the radiation (I < 23 nm in the present case). In metals a 
simple expression for c is given by the Drude formula^: 



e-1- 



(li^ + iojy ' 



(4) 



In this model, a variation in R is ascribed to a change in cOp or y when the photon pulse 
impinges. However, as shown in Fig. 1, in the used co range the dependence on cOp is 
much more significant than on y. So we arbitrarily set Hy =1.5 eV, in analogy with 
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Figure 5 | Ig/ /p^j^-values as a function of Fin the low-Frange for hm = 
18.9 eV. 
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Figure 6 | t- dependence of G{t) (left scale; blue line) and of R{F, t)/R{Q) 
(right scale; black/red and full/dashed lines are Hco = 18.9/20 eV and F = 
20/5 J/cm^, respectively). 

what was found for Al (where y = 0.1 ■ £0^(0)'^) since a reference value for Ti is not 
available in the literature, and we drop a possible P-dependence of)'. We also stress 
that within the accuracy of the present experiment the results are only marginally 
affected by the choice of the /ly- value within the 0.5-2.5 eV range, which essentially 
includes the ?i)'-values observed in a large number of metals^^"^', as also shown in the 
Supplemental Information. For this reason the observed increase in R was related to 
an increase in cOp and then in through Eq. (1). As quoted in the main text, in order 
to provide a quantitative estimate we assumed that is proportional to E and we 
considered both the time (t) and transverse spatial {x, y) profiles of the FEL pulse. 
These result in a (f, x, y)-evolution of the reflectivity {R{F, t, x, y)), which can 
be calculated using Eqs. (3) and (4) assuming cOp{F,t,x,y) — cOp{0){l+AE{Fj,x,y)), 
where 

E{Fj,x,y) / dt' S(x,y)G{t'){l - R{Fj\x,y)) (5) 

^ J — Tj 

is the spatial profile of the energy density at time t deposited in a sample layer 
(involved in the reflectivity process) of thickness d< L below the surface, where G{t) 
and S{x, y) are the normalized time and space intensity distributions. These are 
approximated by Gaussian functions with cr^ = cr^ = 5 fim and = 45 fs. was 

evaluated as Gi — a seed j \/N^^, where (7^^^^ = 70 fs is the time duration of the seed laser 

pulse (measured with an optical cross- correlator) and N — 4 is the harmonic number 
at which the FEL radiators were tuned. Representative i-traces of R{F, t)/R{0) — (1/ 
R{0)) \ dx\ dyS{x,y)R{F, t,x,y) are shown in Fig. 6 along with G{t). It is worth stressing 
that the possibility to approximate G(f) with a well defined function (which ultimately 
permitted us to carry out the present data analysis) is one of the advantages related to 
the use of seeded FELs. Finally, the factor 1 — R{F, t, x,y) in Eq. (5) accounts for the 
reduction in the amount of energy available at later times for further excitation due to 
the f-increase in ojp (and thus in R) within the FEL pulse duration. Finally, the 
observed R{F)IR{Q>) is simply given by the space-time average of R{F, t, x, y), i.e.: 

The mean ion charge (or equivalently the average free electron density) can be 
calculated as well through the definition of cOp (see Eq. (1)), i.e.: 

^^^^ = ^ = dtS(x,y)G(()co^(P,t,^,7). (7) 
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